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I— 1 1 Abstract 

^^ , Effects of spin-orbit coupling (SOC) on metal-Mott insulator transition 

C^ I (MMIT) and spin exchange physics (SEP) of two-component Fermi gases 

in two-dimension half-filling square optical lattices are investigated. In 
the frame of Kotliar and Ruckenstein slave boson and the second order 
perturbation theory, the phase boundary of paramagnetic MMIT and 
r^ ' spin exchange Hamiltonian are calculated. In addition by adopting two 

mean-field ansatzs including antiferromagnetic, ferromagnetic and spiral 
phases, we find that SOC can drive a quantum phase transition from 
antiferromagnet to spiral phase. 
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1 Introduction 



In a crystalline solid spin-orbit coupling (SOC), which occurs naturally in systems with 
broken inversion symmetry and makes the spin degree of freedom respond to its orbital mo- 
tion, is responsible for many interesting phenomena, such as magnetoelectric effect [HEIIH], 
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visionary Datta-Das spin transistor [HIS], topological insulator [6l[7] and superconductivity 
[HllH]. Taking topological superconductivity for example, it has been predicted to occur in 
superconductors with a sizable spin-orbit coupling in the presence of an external magnetic 
field [lOl HTJ [121 El E] • In these systems the transition to topological phases requires that 
critical magnetic field is much larger than the superconductivity gap above which an s-wave 
superconductor is expected to vanish in the absence of SOC. It is SOC that competes with 
a strong magnetic field to give rise to a topological superconducting phase. 

As is known to all that ultracold atom systems can be used to simulate many other 
systems owing to their many controllable advantages and operabilities [151 [161 I17j . Cer- 
tainly the simulations to SOC, which are generally equivalent to produce non-abelian gauge 
potential with optical |18 [ 1191 [20] or radio- frequency fields [211 . are also possible and have 
been realized in a neutral atomic Bose-Einstein condensate (BEC) by dressing two atomic 
spin states with a pair of lasers [22] . Motivated by such a pioneer experiment and a practical 
proposal for generating SOC in ^^K atoms [23], BCS-BEC crossover in the two-component 
Fermi gases with SOC have been widely studied [2l[25l[26l[271[28l[29l[3ni[3l]. 

By contrast in this paper we consider repulsive two-component Fermi gases with SOC 
in a two-dimensional square optical lattice, and are interested in the effects of SOC on 
metal-Mott insulator transition (MMIT) and spin exchange physics (SEP) at half filling. 
Essentially MMIT of two-component Fermi gases without SOC in an optical lattice has 
been realized experimentally [16j and is driven by the competition between hopping term 
and on-site interaction in the frame of one-band Hubbard model. When the hopping term 
dominates, atoms conduct freely in a lattice and the system is a metal. Gradually adjusting 
on-site interaction to an extent that the gain of the kinetic energy cannot offset the increase 
of potential energy, on-site interaction forbids the hopping of atoms and the system evolves 
into Mott insulator (MI). In a MI, SEP is correctly described by quantum antiferromagnetic 
Heisenberg model, which is known from the famous t — J model [32j expected to offer a 
mechanism for high temperature superconductor. Physically this effective antiferromagnetic 
coupling between the nearest-neighbor spins comes from Pauli exclusion principle and the 
fact that the hopping of a particle cannot change its spin. Thus to minimize kinetic energy 
the nearest-neighbor spins must be antiparallel. In the presence of SOC, it has two effects 
on atom hopping. On the one hand SOC can make atoms move from one site to another 
site and corresponds to an effective hopping term, so it definitely has important effects on 
MMIT in view of the above statement. On the other hand, the effective hopping induced 
by SOC is spin-flipped to support the nearest-neighbor spins parallel. From this viewpoint 
SOC also dramatically changes SEP. Furthermore it is likely that when the strength of 
SOC is beyond certain critical value, the system will show a quantum phase transition from 
antiferromagnetic to other magnetic states. 

This paper is organized as follows. In section 2, we firstly derive Hubbard Hamiltonian 
in two-dimensional square optical lattices with SOC, then by using Kotliar-Ruckenstein 
(KR) slave bosons [33j we investigate paramagnetic MMIT, i.e. Brinkman-Rice phase tran- 
sition [34] . In section 3 under the limit of large on-site repulsion and using the second 
order perturbation theory spin exchange Hamiltonian is obtained. By making mean field 
approximations we find that the ground state of the system is either antiferromagnetic or 
spiral depending on the relative magnitude of hopping term and strength of SOC and a 
quantum phase transition happens between them. The conclusions are given in section 4. 



2 Metal-MI Phase Transition with SOC 



The Hamiltonian of the system we consider is 
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where a Fermi atom of mass m for spin a is described by the field operators ^a(r) and 
Vol{^ is optical potential for two-dimensional square lattices. A, g{> 0) and a represent 
the strength of Rashba SOC, two-body contact interaction and Pauli matrix respectively. 
When temperature is very low and filling factor is not too high, all atoms are constrained 
into the lowest band of the optical lattice. Expanding the field operator in terms of the 
Wannier functions ^ai''^ = '^iO.iawir — Ri), where aia is the annihilation operator for 
an atom of spin a in site Ri, and only retaining on-site interaction and nearest neighbor 
hopping, we find 
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where the hopping term t"- is a 2 x 2 matrix and its elements are tj 
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These parameters t, Tf •. 






and U are related to the Wannier 



function as follows 
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and the symmetry of Wannier function, Ff •, F^^- satisfy the 



Prom above expressions ((3]) and the symmetry of Wannier function, F; 
relations F? • = T\ ■ = — F^^ = — F^j. Moreover Ff • = if i,j are nearest neighbor along 
y direction and F^ ■ = if i,j are nearest neighbor along x direction. For convenience the 
parameter F is defined F = |Ff J = |F^ I to represent the strength of SOC. 

It is well known that MMIT is a phenomenon of strong correlation. In terms of strong 
correlation, apart from some numerical methods, such as dynamical mean-field theory [35j, 
a few analytical methods are also available. The first is Gutzwiller variational wave function 
|36j . In this method to make the calculation tractable, one has to introduce the Gutzwiller 
approximation which is basically at the mean-field level. Although this method is successful 
to predict the existence of MMIT, it still has some disadvantages from variational and 
mean- field approximations. Another method is KR slave bosons ^33j. It exactly reproduces 
the results of Gutzwiller approximation at the saddle-point level and can be improved 
systematically by considering fluctuations around the saddle point [37]. Hence below we 
adopt slave bosons to study MMIT with SOC, although our results is at the saddle-point 
level. 



For two-component Fermi gases, the Hilbert space for every lattice site i consists of four 
states |0 >j, \a >«= a|„|0 >j and | tii>«= Q^I-i-oLlO >«• ^^ the representation of KR slave 
bosons, in addition to original fermions, a set of four bosons e, d, Pa for every lattice site 
are introduced so that |0 >i= e\\vac >, \a >i= PiaoJal^'^'^ ^« ^^"^ I tii>i= dia-^a-Avac >i, 
where \vac > is the vacuum state after introducing slave bosons. It is easily found that 
ejej, dldi and pJaPia represent the projectors on the empty, doubly occupied and singly 
occupied site. Due to the fact that the introduction of bosons enlarges the Hilbert space 
of every site to contain some unphysical states, such as e\e\\vac >i etc., we must impose 



three constraints ejci + pl^Pia + d\di 



1, at a 



and considering above constraints the Hamiltonian ([2]) is reformulated into 
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As claimed by KR, the substitution Zia for zjq, ensures zJ^Zjis = 1 to recover the results in 
the limit U = at the saddle-point approximation. 

The partition function Z can be written as a functional integral over the fermion and 
boson operators 
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where the Lagrangian C{t) is 
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and //, XI , XlJ are the chemical potential and Lagrange multipliers, respectively. 

Assuming uniform and static boson operators and Lagrange multipliers, i.e. at the saddle 
point, one can integrate over fermion operators and obtain for thermodynamic potential of 
single site 



n = A«(e2 + ^2 +pI- 1) + Ud^ - Ai2)(d2 +p2) _ J_^in[l + e-/5^-] 
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with E, 



ka 



Cfct + e^ + av/(efct " ^h)^ + ^^f^f^l /2> ^fca = ekzl-fi+X^a', €k = -2t(cos k^a+ 



cos kya), Ffc = 2rw sin^ /c^^a -|- sin^ fcyO. It is to be noted that a, N are lattice length and 
the number of lattice site, and wavevector k belongs to two dimensional Brillouin zone. At 



this time the seven parameters e, pa, d, X^^' and \a are obtained by minimizing Q, and 
the chemical potential at half filling by thermodynamic relation — gr^ = 1- These equations 
are called saddle-point and number equations. 
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Substituting above relations into saddle-point and number equations, one still has two 
equations satisfied by fj, and d 
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At zero temperature, in the frame of KR slave bosons, d^ = corresponds to the 
vanishing of the number of doubly occupied sites and indicates that the system is undergoing 
a MMIT. From this criterion one has numerically solved the equations (|10p . The numerical 
results suggest the chemical potential is still fixed at /x = U /2. Hence at zero temperature 



^ = ^ E {(^fc + rfc)e[-(efc + Tk)] + {ek - Tk)Q[-{ek - Ffc)]} . 



(11) 



and the phase boundary of MMIT is 



u = -se. 



(12) 



where 0(x) is Heaviside step function. Without SOC, ^ = -^ X^fc £fc0[— efc] and the phase 
boundary (|12p is the same as the result in |33j. In Fig.l the phase boundary of MMIT 
is shown. From Fig.l, very explicitly SOC stabilizes the MI, which is consistent with the 
fact that SOC can be regarded as an effective hopping term. Besides instead of adjusting 
t MMIT can also be driven by changing SOC, so one has found another way to realize the 
MMIT. 



3 Spin Exchange and Magnetic Phase Transition with SOC 

As demonstrated in the section 2, at half filling when U » t and [/ >> F, the hopping of 
atoms are forbidden and the system evolves into MI with spin S = ^ for every lattice site. 
In the MI we could regard t and F as perturbations. In the limit of t = F = the energy of 
the system does not depend on the spin orientations on different sites. When t, F are finite 
but small, we expect that we still have spin 5 = 2 in each site, but atom hopping processes 
induce effective interactions between these spins, usually called spin exchange interaction 
|,38j . To construct an effective spin exchange Hamiltonian for this system, we note that 



in the second order in t, T it can be written as a sum of interaction terms for all nearest 
neighbor sites. These pairwise interactions can be found by solving a two-site problem in 
the second order in t, T. 

The ground state manifold for two-site problem with one atom in each site composes of 
four degenerate zero-energy states 

|1 >= I t>i I t>i, |2 >= I t>i I i>j, 

|3 >= I l>i I t>i, |4 >= I l>i I l>j, (13) 

with i,j labelling two sites. The first order perturbation theory takes us out of the ground 
state manifold and can be neglected. In the second order atom hoppings can connect all 
four states by two intermediate states |5 >= | tii>j |0 >j and |6 >= |0 >i \ '[,],> j. To 
find spin exchange Hamiltonian we need calculate all matrix elements |i38j 

where states \a >,|6 > and \c > respectively belong to ground state manifold and in- 
termediate states with E^ representing eigenenergy of corresponding state in the zeroth 

order. The calculation is very direct and when i,j are nearest neighbor along x direc- 

2r^ pa; 2tr^ 

tion, we have Mn = Mu = M41 = M44 = '^ '% M12 = M21 = M24 = M42 = -^, 

Mi3 = M31 = M34 = M43 = -jj^, M22 = -M23 = -M32 = M33 = -^. Accord- 
ing to spectral representation of an operator, magnetic Hamiltonian of two-site problem is 
Hi J = J2a,b \^ > ^ab < b\- Making substitutions |1 >— >• aLaL, |2 >— >• a-^a-,, |3 >— )■ oL oL, 
|4 >— )• oL oL and using algebra of spin operator Si = ^ala^apiip^ we get 
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By the same procedure, when i,j are nearest neighbor along y direction we get 
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Thus spin exchange Hamiltonian of the whole system is 

Hse = E ^*J- (17) 

Some comments about ()17p is following. If F = Hge describes isotropic quantum 
antiferromagnet, consistent with t — J model. When F 7^ 0, main effect of SOC is to break 
spin conservation by two ways, one of which, corresponding to the second term in (jlSp and 
()16p . flips one spin of two nearest neighbor sites, while the other flips simultaneously two 
spins corresponding to the third term in (J15p and (|16p . Thus the antiferromagnetic state 
will be unstable when the strength of SOC F is beyond certain critical value. 

Now we decide the ground state of the system at mean-field level. This corresponds to 
regard quantum spin operator Si as a classical vector. The first mean-field ansatz including 
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ferromagnetic and antiferromagnetic states is that spin configurations in two sublattices of 
a square lattice take different values specified respectively by coordinate angle {-&, ip) and 
(7, 5) with < t?, 7 < TT and < ip,6 < In, the mean-field energy scaled by U is 



E = 8N^ 



(t — r )cos'/?cos7 + t sint?sin7cos(93 — (5) , (18) 



where a 2N x 2N lattice is assumed and t = t/U, T = T/U. Easily found that energy only 
depends on the difference of if and 6, for convenience we can choose 5 = 0. Owing to factor 
sin 1? sin 7 > 0, the minimization of energy leads to if = it. Minimizing energy about t?,7, 
we get 

(P — r^) sin -d cos 7 + P cos ■!? sin 7 = 0, 

(? -f^)cosT?sin7 + ?sim9cos7 = 0. (19) 

Equations (|19p have two sets of solution ?? = 7 = and t? = 7 = tt/2. The first solution 
corresponds to ferromagnet along z direction with Efe = 8N'^(lr — F^) and the second 
corresponds to antiferromagnet along x direction with Eaf = —SN'^lr. If EpE < Eaf 
ground state is ferromagnetic, on the contrary ground state is antiferromagnetic. Thus this 
mean-field ansatz predicts a phase transition from antiferromagnet to ferromagnet and the 
critical point is EpE = Eaf, i-e. T = v2t. 

The motivation of the second mean-field ansatz comes from i = limit in spin exchange 
Hamiltonian (117p . Letting t = the classical spin configuration minimizing energy satisfies 
three conditions: (1) z components of all spins are equal; (2) for a random chain along x 
direction x components of all spins are equal but y component must be alternating; (3) for 
a random chain along y direction y components of all spins are equal but x component must 
be alternating. Such spin configuration, which we call spiral phase and shown in Fig. 2, is 
permissible in a square lattice. From above three conditions if coordinate angle (0,0) of a 
spin in the lattice is specified, energy of the system is 

Esp = 8N^(P cos^ 0-f^), (20) 

and its minimization gives rise to 9 = tt/2, Esp = — SA^^F^. Comparing Esp with Epp 
we find that the ferromagnetic state is always a metastable state. As a result phase tran- 
sition predicted by the first mean- field ansatz does not exist, we get a phase transition 
from antiferromagnet to spiral phase with critical point T = t. Fig.l also shows magnetic 
phase diagram in terms of such two mean-field ansatzs. Physically the metastability of 
ferromagnetic state is attributed to the fact that SOC breaks spin conservation. 

4 Conclusions 

In conclusion we have discussed MMIT and SEP of two-component Fermi gases with SOC 
in two-dimensional half-filling square optical lattices in the frame of KR slave bosons and 
second-order perturbation theory. Comparing with the case without SOC, SOC not only 
enlarges the region of MI in the phase diagram and introduces another way to realize MMIT, 
but also dramatically affects SEP due to SOC breaking spin conservation. Importantly 
by adopting two mean-field ansatzs we find that SOC can drive a phase transition from 
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antiferromagnet to spiral phase. Experimentally this phase transition can be observed by 
either adjusting optical lattices to suppress the hopping term or decreasing the strength of 

soc. 
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Figure 1: The phase diagram of a repulsively interacting two-component Fermi gas with 
spin-orbit coupling in a square optical lattice. The solid line is the phase boundary of 
metal-Mott insulator transition, while the dashed line is one of antiferromagnetic-spiral 
phases. 
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Figure 2: Spin configuration of spiral phase is shown with small squares representing lattice 
sites. {9, (p), {9, —4>), {9, it — 4>) and (9, tt + 4>) are coordinate angles of spins. 
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